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We calculate the finite size scaling of the sample-to-sample fluctuations of the free energy of 
the m component vector spin glass in the large-m limit. This is accomplished using a variant of the 
interpolating Hamiltonian technique which is used to establish a connection between the free energy 
fluctuations and bond chaos. The calculation of bond chaos then shows that the scaling of the free 
energy fluctuations with system size A'' is AF ~ A'^'' with \ < n < j^, and very likely M = | exactly. 

I. INTRODUCTION 

Spin glass physics^ continue to be a very exciting and difficult topic. One of tlie important ongoing issues are 
the finite size corrections to thermodynamic quantitiea^rii. Such finite size corrections are usually impossible to 
calculate witfiin replica tlieory (which is otherwise extremely successful for the spin glass and other problems), due to 
the massless modes which are often encountered and which prevent going beyond zero-loop order in a perturbation 
expansion. Nevertheless, many results are by now established for the Sherrington-Kirkpatrick spin glas&i^, either 
numerically or analytically (see references cited above). 

An observable which is particularly interesting is the finite size scaling of the sample-to-sample fluctuations of 
the free energy. This quantity provides a link between two apparently distant fields, namely spin glass physics 
and extreme value statistics^'^. In extreme value statistics, the question is the probability distribution of extremal 
events such as the maximum (or minimum) of a set of random numbers. The classic results of extreme value theory 
state that such extremal events follow one of three possible limiting distributions (the WeibuU, Gumbel or Frechet 
distribution). Recently, a fourth universality class was found, the Tracy- Widom distribution for the smallest (or 
largest) eigenvalue of a Gaussian random matrix^**. Similarly, one could ask what the distribution of the ground 
state energy of a (Sherrington-Kirkpatrick or other) spin glass is, which is a question of extreme value statistics. In a 
statistical mechanics setting, however, it can be generalized to the question of what the distribution of the free energy 
is at finite temperature. This is a very difficult question indeed (both for finite and zero temperature). To keep it 
simpler, we merely ask what the width of the distribution is (i.e. the sample-to-sample fiuctuations) , and this width 
Ai^ will scale in some way with the system size iV, i.e. A_F ~ with an exponent /z. These free energy fiuctuations 
have been considered for the Sherrington-Kirkpatrick model by numerical investigations, see e.g. Refs/llUS There 
also exist heuristic arguments for /i = and /i = I9,ii,i5- _i7^ g^j^^j ^^^q limit n < \ has been shown -^-li^. All these 
results show that the Sherrington-Kirkpatrick model does not fall in any of the four established universality classes 
of extreme value statistics. For a different famous replica symmetric spin glass model, the spherical spin glasses, 
the situation is different. Since its groundstate energy is the smallest eigenvalue of a Gaussian random matrix, this 
model falls into the Tracy- Widom universality class^. This implies that the fluctuations of the free energy scale as 
AF N^/^, which has been confirmed recently by a replica calculationii. 

In this work we consider the fully connected m component vector spin glass in the limit of large m. This model is 
known to be replica symmetrio^^, which suggests that its free energy fluctuations might fall into the same universality 
class as the spherical model. Furthermore it might be expected that due to replica symmetry the fluctuations in 
this model are simpler to calculate than in the Ising spin glass. Unfortunately, this hope is not entirely justifled, as 
we will see. As mentioned above, it is usually impossible to calculate subextensive quantities within replica theory. 
We circumvent this problem by using a connection between the sample-to-sample fluctuations of the free energy and 
bond chaos, which was derived for the Ising spin glass^^, and which we generalize here to vector spin glasses. This 
connection allows us to calculate the sample-to-sample fluctuations by calculating bond chaos instead, and this is 
possible within a large deviation approximation, combining techniques from Ref. ^23! and [2^ . As it will turn out, 
the large deviation approximation is not good enough to obtain the final answer, and we must resort to additional 
resources, such as our previous knowledge of the scaling properties of the large-m model^^"— and some additional 
scaling assumptions. To forestall our main result, we obtain < ^ < ^ and most likely M = f exactly. This shows 
that the large-m spin glass does not belong to any of the four universality classes and probably to a different class 
than the Ising spin glass (unless it so happens that e.g. M = i for both the Ising and the large-m spin glass, which is 
not entirely ruled out by our results but which we deem unlikely) . The result l^, — \ was first suggested in Ref. [25l . 
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However, the derivation required the existence of a gap in the eigenvalue spectrum of the inverse susceptibihty matrix 
which was later shown not to exist^^. Here, we provide a different explanation for this value of /i. 

This paper is organized as follows. In section|ll]we will show how to derive the exact connection between the sample- 
to-sample fluctuations of the free energy, A.Fn, and bond chaos for the m-component spin glass using interpolating 
Hamiltonians. We will see that the fluctuations only depend on the finite size scaling of averages of powers of the 
link overlap ql between two copies of a spin glass with different (but correlated) disorder. The averages of the link 
overlaps will be calculated in section Hill bv calculating the probability distribution of the link overlap, Pe{qL), where 
the parameter e measures the degree of correlation. We will use these results in section HVl to derive the exponent ji 
of the free energy fluctuations. 



II. CONNECTION BETWEEN THE SAMPLE-TO-SAMPLE FLUCTUATIONS AND BOND CHAOS 

We use the technique of interpolating Hamiltonians^^ in order to derive an exact connection between the sample- 
to-sample- fluctuations of the free energy, AFn, and bond chaos in the m-component spin glass. To do so, we adapt 
the method from Ref. which was developed for the Ising spin glass, to the m-component vector spin glass. The 
Hamiltonian of the m-component spin glass is 

i<j 

where Si are m-component spins and are independent Gaussian random variables with unit variance. The spins 
are assumed to be normalized in such a way that sf — m. 

The main idea for calculating the sample-to-sample fluctuations of the free energy is to use the following two 
Hamiltonians, 



'Ht^ -\j^^-—^^J^JS,Sj- J^^J[jS,Sj and 



N A^"'^"'"-' V N ■ 

t<j i<j 



1 -T 

i<j ' i<j 



with < i, T < 1 and J[j and J^j additional independent Gaussian random variables with unit variance, and express 
the free energy fluctuations in terms of them. These Hamiltonians interpolate between a spin glass with a given 
set of coupling constants {Jij} for t = (or t = 0) and an identical spin glass with a different, independent set of 
coupling constants {J'ij} or {Jij} for t = 1 (resp. t — 1). For all < i, r < 1 we have the same type of spin glass but 
with coupling constants — tJij -f \/tJ[j (or \J\ — rJij + y'rJ^ ), which are still Gaussian random variables with 
unit variance. We denote the disorder average with respect to all coupling constants {Jij}, {J'ij} and {J'ij} as E ■ ■ ■ . 
Thermal averages will be denoted as (•••). 

The free energy fluctuations can be written as /3^AF^ = i?(log Zi — log Zo)(log Z[ — log Zq), with Zt and Z'^ being 
the partition functions with respect to the Hamiltonians Ht and H'^.. This can be seen by writing logZ = —f3F and 

using the independence of the different sets of bonds for EFiFq = EFqF[ = EFiF[ = F and EFqFq — F^ (the 
overbar denotes a disorder averaged quantity of the original system of interest, i.e. Eq. JT])). We use the idea in Ref. 
to represent this expression by integrals in the form 



(3) 



To calculate the right hand side, we follow Ref. Il9| and obtain 
;±llogZt^logZ' = ^ 



E— log Zt— log Z'^ = TfiA72^E (((«^^j)(^fcsO>t - {s^Sj)t{SkSl)t){{{s^Sj){skSl))r - {s,Sj)r{SkSl)r) 



ijkl 

13' 



(i?E {{s^SJ)t{s^S,)r) - n) . (4) 



The subscript t or r on the thermal averages indicates whether the average is to be taken in a system with Hamiltonian 
-Ht or n'^. 
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We introduce the link overlap between two replicas with potentially difTerent coupling constants as q]^ = 
N{N-i) J2^<J{siSj)t{siSj)r = jv(jV-i) ( ( ( E» g^;!^ V) " m^Af) . As wc will shortly need not only 2 but up 
to 4 replicas, we label the replicas with upper indices (1) to (4), where replicas 1 and 2 have Hamiltonian Ht and 
replicas 3 and 4 have Hamiltonian H'^. Analogously to q]^, we define the other possible link overlaps q}^, qf^ and qj^. 
Then equation Q can be expressed in the following way 

The notation [. . . ] in the last line stands for the average taken with the bond averaged probability distribution P^{qL) 
of finding a given link overlap qL- The parameter e indicates the statistical "distance" between the sets of bonds in 
the two replicas involved and will be defined in detail below. In principle, one would have to consider simultaneous 
multi-replica overlaps (such as, for instance, the term E{qj^q]^)). Fortunately, it follows from replica symmetry, 
which holds for the m-component spin glass in the large-TO limit, that the corresponding joint probability distribution 
P}^^ {qj^ , qj^) factorizes into P}^'^{q]^, qj^) = Pe{q]^)Peiqj^)- A similar statement holds for the four-replica probability 

distribution P,1234(^13^^24) 

Up to now the overlaps ql'' depend on the two parameters t and r. However, the only important quantity is how 
much the sets of bonds in the two replicas a and b differ. For example, when i = t = 0, the bonds in replicas 1 
and 3 (say) are identical. On the other hand, when t = and t = 1 (or vice versa), the bonds in replicas 1 and 3 
are completely uncorrelated. The degree of correlation between the two sets of bonds can be measured by the one 
parameter e defined by ^^^^^ = y/1 — t\J\ — t, as shown in Ref. 19. When e = 0, the bonds are identical; when 

e = CXI, the bonds are completely uncorrelated. 

The parameter e can be used to eliminate the integration variable r from Eq. by a variable transformation. The 
integral over t can then be evaluated exactly, and as a result, the connection between the sample-to-sample-fluctuations 
of the free energy and bond chaos in the m-component spin glass is found to be (see Ref. 19; for details) 

/J^AF^ = (iV_^ &/2(.)([(<zF)^] - \qff) + li^L-lM! deg,{eM% 

- ^21 + ^ 122 (6) 

with the nonnegative functions /2(e) — '^'^^(^^^2^2 ^ and (72(e) = ^^^^^^^^/t- Note that bond chaos enters Eq. © 
through the measure of distance e of the bonds of the two replicas between which the link overlap q]^ is calculated. 
The integrals I21 and I22 will be calculated below. 

The analog of Eq. 1^ was called "second route to chaos" in Ref. Il9j, hence the first index is 2 on I21 and l22- In 
addition to this result, however, there is a another exact relation between the fiuctuations and bond chaos (the "first 
route to chaos"). It stems from using only the first interpolating Hamiltonian of equation ([2]), "Ht, and the relation 
20^lS.Ffj — E{\ogZi — logZp)^, which is easy to derive. Proceeding similarly to above, one can prove the following 
equality: 

/3^AF^ = -^^^^ f deMemqir] kF?) + i^^^^ deg,{e)[qf] 

^11 + 1 h2. (7) 

The only difference between this equation and Eq. ([5]) is the minus sign in front of the first term and the weight 
functions /i(e) = (j^^2y arcsin ^-^^^2 ^nd (?i(e) = ^]^_j_J^3/2 arcsin ^-^^^2 in the integrals In and I12 instead of /2(e) 
and 32 (e) as in /21 and /22. The minus sign of the first term implies that the second term is an upper bound of 



III. CALCULATING P,XqL) 



Eq. ([SJ involves moments of the link overlap taken with the probability density P^iqh)- This function can in 
principle be calculated by taking two replicas with bond realizations drawn with parameter e and constraining the 
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replicas to have link overlap ql- This constrained system has free energy F^^j^ql), and the (non disorder averaged) 
probability density Pe,j{qL) then follows to be 

Jo "9l exp(-pi*e,j((7ijj 

Finally, P^j{qL) must be averaged over the disorder. 

Unfortunately, this task is too difficult in general. Instead, we will calculate only the disorder averaged extensive 
part of the free energy, denoted by Nmf^{qL), and the probability density defined by 

po/ X ^ exp(-/3iVm/,(gL)) , . 

'^'^^^ dq'^eM-^NmUq'^))- ^> 

This is the large deviation approximation to P^lq^). Averages taken with respect to P^{qL) will be denoted by [. . . ]o. 
The tail of this distribution will be the same as that of Pe{qL), but in general there will be deviations. In fact, this 
is the point where this paper differs most from Ref. Il9| since in that publication, the difference between the true and 
the approximative distribution was only quantitative, whereas here it is substantial. We know that for e = the link 
overlap distribution consists of a 5 peak at the Edwards- Anderson value and elsewhere since the large-m spinglass 
is replica symmetric. We will see below, however, that we do not observe this peak in P^(cii,) at all. It must therefore 
be generated from the finite size corrections to the extensive part of the free energy. 

Hence the finite size corrections are very important in this calculation, but we have no direct way of calculating 
them. In order to overcome this problem, we will show that P^{qL) will be valid for e larger than some crossover 
value, and we will use additional arguments and simulation results to fill the gap for smaller e. 



A. Replica calculation 

The first task is to calculate the extensive part of the disorder averaged free energy Nmf^{qL) of two replicas 
constrained to have link overlap q^- To this end we calculate the partition function Z^^j^q^) of this two replica 
system. This system has vector spins with (s^)^ — m for every spin i = 1,...,A^ and the two real replicas 
X — (10), (2e) (this notation denotes the replica number in its first entry and the value of e in its second entry). These 
two replicas differ in their coupling constants KiAf) — . ^ K^, H , ^ K', , where if? and K', are independent 

Gaussian random numbers with unit variance, by choosing ^ = for the first replica and ^ = e for the second replica. 
Furthermore, the two replicas have link overlap g^, which is enforced by a (5 function in the partition function as 
follows: 

zM = Tr.-(.{,. - j: ( x: - -^)) 

X -P (;4 E Kmr^^'' + ;4 E (10) 

We follow the replica calculation of Viana^"^ for the n times replicated partition function, write the ^-function in an 
integral representation with parameters Zq,, and introduce the traceless tensor T'^'^^s -j-q separate the diagonal part of 
QZ (with Qa^ = 2 ^ ;<_f + 2 Ea,M<. T^^" + Ea,^ (Q- + ^a"") )• After evaluating the trace and regarding 
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only terms up to third order in the tensors Q, T and R we get 



m + 2' 



TO + 2 



3! ( E Qa/sQ p'^Q'ya 



(to + 2)(to + 4) 



ErpjAl/rpUprppfj. 



(m + 2)(m + 4) 



E(^r)' 



3m 
TO + 2 

3771 
TO + 2 



E Qr.Q^^.T^r + 3 E Qr,^7.^r. + 3 E q^x^t. 



a</3 



(A1<>^),P 



with T = i(l 



M ^' - 1 



(1 



- 1 , r' = 2to — 1 , o; = 1 and the matrix A'^'i is given by 



A = 



(11) 



(12) 



with the above separation for and only T^^ on the diagonal of the matrices Q. We split the tensor R into its 
diagonal matrices p'^" — R^'^ and the rest. This equation is solved by a saddle point integration for the tensors Q'^'j^, 

, 1^ n -4- In r-\ i^ -r* -rvi /^-i- y^v^fr ^ 



Tj^" and R'^^q and the parameters Zo 



B. Solving the saddle point equations 

Solving the saddle point equations is the remaining task to derive the characteristic form of the overlap distribution 
Pe{qL)- As the TO-component spin glass was shown to be replica symmetric in the limit to od^, we calculate the 
saddle points in the replica symmetric case. The ansatz for a replica symmetric scenario is as follows^^i^: 

= P^iii^^ pT = PdS^iv and = z (13) 

To briefly justify these equations, one has the usual interpretation Q^^'^'^' = (sq^^'s^^j^^) from the replica calculation. 
Due to the isotropy of the model, averaged quantities like this reduce to x5^j_i, with some mean value x, depending on 
the quantity at hand. Thus EZ^j{q) is given by 



'"j(9l)0C J (^^-I)"dze^27Vn9^^-r„^«.y rfA^- gxp (iV [rTO77(77 - 1)0^ + /^^(n - l)p2 + !^p2 

1 /2z 1 ^(-1) 2 , o ^"("-l)("--2) „3 

- "-^(^ + 7f^) + 2 

+ ujn{n - 1)(77 - 2)jnQP^ + 2ujn{n - l)mQPpd ^ (14) 

This leads to four different saddle point equations (we use Nm as the large parameter) 

= tQ-Q^-P^ + Ppd (15) 
= T'P-2QP + Qpd (16) 

0=p,-^p,C^ + ^J=Y''^-2QP (17) 
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where line 3 and line 4 can be combined to 

Pd - 2QP = ./^. (19) 
V m 

In the following, we substitute — ?> q. 

From equation the free energy is (after taking the usual replica limit n — ?> 0) 

/?/,((?) = const. - ^ ~ + rQ2 + /P^ _ _ 2Qp2 ^ 2QPp, (20) 

1. Above and at the crittcal temperature 

Above (r < 0) and at the critical temperature (r — 0), the Ising spin glass is replica symmetric, just as the m- 
component spin glass. Therefore the solutions of the saddle point equations are the same as in Ref. [2^. Inserting 
them into the free energy (equation (|20l) ) gives 

above Te: /3/,(g) = ^ (l - — £L=) + 0(^4) (2i) 



atT,: \ ^ \1 ' : (22) 



At Tc the probability distribution Pe(g) oc e consists of two parts with different dominating exponents 

depending on e being smaller or larger than N^^^^: 



P^(.l)^< „ (23) 



'Nm3^ iV-1/6 ^ e 



2. Below the cntical temperature 



The equations ([T5l). (116)) and (1191) can not be solved for general q and e. Due to this, we calculate corrections to 
the two following solvable cases, q = and e = 0, perturbatively in various limits. The exact result for g = is 

P = 

Q - r (24) 



To find the solution for e = 0, we rewrite equations (|15p and (|16p in terms of the new variables a = Q + P and 
b = Q-P 

{T + pd)a-^^ia-b)-a^ ^0 (25) 

{T-pd)b-^^{a~b)-b^ = (26) 

with the solution (e = 0, so r = r') 

Q = r (27) 

P = P. = ^, (28) 

for a 7^ and 6 7^ 0. With equation (P(I| we get /3/o(9) = /3/o = l''''^- We will use this free energy as the reference free 
energy as it is the energy of the unconstrained system. For q > t ~ 2t^ — qEA the solution with 6 = maximizes the 
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free energy (the solution with a = is an unphysical one). In terms of Aq = q — [t— 2t^) it is (to lowest order in Ag) 



= P = 



2 



Vd = T+-^, (29) 

and the difference to /3/o is 

A/o(<z) = coAg^ (30) 

with Co = 6(i-2t)'' • -l^ow we will calculate corrections to the first solution, in various different limits of e and q. 



In the limit e — >■ cxd, i.e. r' — >■ — c» Eq. (ITCl) yields the solution (to leading order) 



P=-^+0(^) (31) 



For equation (fT?)) this leads to (again in leading order) 



2 

0=(t + ^)Q-Q' (32) 
r 



For r' = —oo,pd is equal to g and the free energy difference to /o is 



2 2 

/J/oo(g) - /3/o = qPd-^+ Oip'd) - |- + 0(g4) (33) 



4- Perturhative solution for <^ q <^ 1 



We use the saddle point equations for a ([25]) and b (|26p , insert a — + ai, b = bo + bi and pd ~ pdo + Pdi , where 
the index denotes the undisturbed solution (e — 0) and the index 1 the first order corrections to it. The results are 
{At = t - r' > 0) 

"1" "M^a - 4T2(l-2r) 

2 

and the correction to pd is Pdi = For the free energy these results yield: 

/?/.(?) - /3/o = ^^^'9' + (^(e'?') (35) 

LOT 

5. Perturhative solution for g ^ <C 1 

In this case the suitable reference solution is the one with q = 0, equations (|24l) . We again introduce corrections 
(Q — > Q + AQ) to this solution and combine equations (fTH]) and (|19p in lowest order to 

r'AP- 2QAP + Q(Ag + 2QAP) = 0, (36) 

where Ag is understood to be the correction to g = 0. The solution is 

T'-2T + 2r2 l-2r 4r(l-2T)' ^ ^ 
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AQ = and 



Ap.-A V_2r + 2r^ ^ T^^ " 2(T3^- ^^^^ 
The result for the free energy then is 

^feiq) - /3/o = „ . ^V + (39) 

lDr(l — 2r) 

The restriction <C 1 is basically unnecessary and we could calculate the free energy in the limit q ^ min(l,e^). 
However, since the regime of large e will not contribute to the sample-to-sample-fluctuations we neglect it already at 
this point, and write (3fe{q) — Pfo{q) — f{e)q^ with /(e) = iqt(i-2t) small e. 



C. P0(g) 

Although we are interested in t he pro bability density of the link overlap qL, the more practical quantity for the 
calculation turned out to be q = \/'qjJm. We therefore formulate our results in terms of q instead of qL for the time 
being. This is not a serious restriction since we have the simple relation 

[(9")"]o = m" / dqq'^^P^iq) m"[q'-]o. (40) 



The probability distribution P^{q) ^ e~^"^0ifei<i)-fo{q)) divides into 4 parts, depending on the range of e. For 
e ^ N~^/^ the contribution to the probability distribution of g of both equations (1351) and (|39l) (both with /e(q) ~ ^.^q"^) 
is negligible, therefore it is approximately a constant in that range for all q e [0, qsA]- Eq. (|30|) implies that P^{q) has 
an exponentially decaying tail for q > qEA with e-^"^<^o(q-qEA) define a function Q{q — qsA) which combines 

both properties. Instead of this plateau in P^{q) there should be a 5 peak at q — qEA which we do not see in 
our calculation. This is due to the fact that we have neglegted finite size corrections to the free energy which are 
dominating in this regime and which we will implement in the next subsection. 

The two solutions we found perturbatively in Eqs. (I35|) and (j39p both produce a probability distribution of the 
form g-^"i0'^^<^ 9" with different constants Cx, but hold in different ranges of e, depending on the relation of and 
q. The order of e determining the crossover from one regime to the other is where is of the same order as q and 
Ne^q^ (in the range iV~^/^ ^ e ^ eo) is of order one such as to be the dominating part of the free energy. This leads 
to e ^ N^^/^ as the crossover value. Thus we get the final result 

\Q{q~qEA) e^N^^'^ 

1 e-^mmg- < 



with ci = and C2 = ,} „ ^ . 

^ 16r ^ 16t(1 — 2t) 



D. The overlap distribution at e = 

In order to account for the missing S peak in P^{q) for small e, we have to consider the finite size corrections to 
the free energy, which are impossible to calculate. However, not all is lost since at least we know that at e = they 
grow as N^~'^ with y = 2/528,^ and we expect that this scaling is independent of the value of e. We can therefore trust 
our result derived above when NmP{f^{qEA) — fo) ^ N"^^^- Since qEA = C>{1), it follows that this is the case when 

Ar3/5 ^ ^-1/10 ^ ^l^gj^ J^J^yg 

Upsiq-qEA) e«7V-i/io 

Pe(g) oc<^ e-^^^^*9' iV-i/io^g^£g (42) 

[e-^/W9' eo<e. 
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FIG. 1: Finite size scaling plot of the variance of the link overlap at temperature T — 0.6. See text for explanation. 

The function SpsiQ ~ Qea) stands for a function which goes to a (5 peak in the thermodynamic hmit. Note that the 
new regime completely replaces the first two regimes we calculated in Eq. (|4ip . 

What we need in the actual calculation of the sample-to-sample fluctuations of the free energy are the expressions 
[(^F)^] ~ ^-iid [qj^]- In the regimes where we have an explicit expression for P^{q), we can calculate this directly 

(see below). In the regime just found, however, we do not have this information available. We must therefore resort 
to other methods and use a finite size scaling ansatz of the form 



(43) 



for e <C iV ^^^^ and with a scaling function J-„i{x) whose properties will be discussed below. Similarly, we assume 
that in this regime [q]^] can be written as 



m 



(44) 



with an exponent p and a scaling function Grnix). These scaling functions and exponents will enter the calculation of 
the fluctuations below. 

We assume to — >■ oo such that the system is replica symmetric. Therefore, the left hand side of Eq. P5)) . which is 
the variance of q]^ /m, goes to for — >■ cc since it is just the squared width of the peak in P^{qL). The exponents 
a and /3 are unknown, but we can obtain a from a simulation at e = by measuring the width of the peak in the 
distribution of qL- This is shown in Fig. [TJ There is, however, a complication involved. In a simulation, m must be 
finite, and for finite to, the thermodynamic limit is a replica symmetry broken phase and the variance o f qV^ w ill not 
tend to 0, thus apparently making a = 0. In order to overcome this problem, we recall that in Refs. I25ll26l it was 
shown that at T = there is a critical number of spin components tiq ~' (with an exponent /i' = 2/5) above which 
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the system does not depend on the number of components any more and is thus identical to the repUca symmetric 
m = oo Umit. We conjecture that something similar happens at finite temperature, i.e. when m ^ , the system 
is in a replica symmetric phase, whereas for m <C iV^ it is in a different phase. For this reason we have plotted the 
variance {Aqj^y of q]^ against x — Nm~^^f^ , such that we can expect replica symmetric behavior for small x and a 
crossover to a constant variance for large x, and this is precisely what can be observed (although the crossover is so 
slow that we do not see the expected plateau yet). For the determination of a only the data for small x are relevant. 
For the simulation, we have implemented a parallel tempering Monte Carlo algorithm using system sizes from 
= 64 to iV = 216 with up to 29 different temperatures in the range of [0.6 : 1.6]. In that way we produced for 
two different replicas (with the same set of coupling constants since we are considering e = 0) at least 64 statistically 
independent sets of spin configuration for every temperature. From these, we calculated the variance of the link 
overlap {Aqj^)"^ — [(g^)^] — ['Z/,'^]^- According to the scaling ansatz above, m~^{Aqj^)'^ ~ {Nm~^^'^ )~"m~"/'' at 
e = 0, such that m°'^f''~'^{Aql^y ^ (TVm-^/^')-". The data in Fig. [T]shows on the one hand that a/^i' - 2 w 1 and 
on the other hand that a sa 1.4. The observation a//x' — 2 « 1 implies a « |. Together with the second, more direct, 
observation of a, the data shows that a > |, and as we will see below, this is all we need to know. 



IV. CALCULATING [g'> AND THE SAMPLE-TO-SAMPLE FLUCTUATIONS OF THE FREE ENERGY 

The remaining task is to calculate [q"']o and insert it into equation ([5]). To do so, we use steepest descent methods 
to write for the regime e with iV-Vio ^ £ ^ (^ith P^{q) = ^6"^^""-^^^^^ and the normahzation constant Afq) 

1 z*^ 

n^ir dgg-e-^^-^V+oC,^)^ (45) 

Jo 



rA 



neglect the term of order q in the exponent (note that we set the upper bound from 1 to oo, which introduces only 

r(^) s 



T( — ) 

exponentially small errors), and get (with Mq — . ^ ^ ) 



1 1'°° 

^ (Wr:^r(!i±i)^ ,40, 

We then have 

With this, we calculate the sample-to-sample-fluctuations through Eqs. ^ and ([7]) by taking the leading order in N of 
every integral into account. The first integrals, J21, (with /2(e) = 2^^ + 0{€^)) separates into three integration intervals 
corresponding to our three regimes. We neglect the range of e > cq, because it gives a contribution of order 1 which 
we are not interested in. For the part e ^ _/V~i/io we have the scaling ansatz for [(gp)^] — [q'^']'^ = rri^ N~°' F„i{N^ e). 
This ansatz must match m^([(7'*]o — [^^lo) from the neighboring regime at the crossover point e — 7V^^/^°, which there 
goes as 7V~^e~^ ~ jsf-^/^ (ggg gqg^ (j47)) and (|48p ). This implies that the scaling function J^m{x) decays as x'^'^ jrr? 
for X ^ 00 with an exponent 7 obeying 7(/3— j^) — | — a. 
This leads to 

^r-l/lO 



hi = / 2e37V-"J-„(iV^e) + / 

Jo JN-^n° 



de 2e^ 



/.jV-i/i" .£0 -5 

^2N-°'-^'^ dxFm{x)x^ + N-^ I de— + ... 

~ const.iiV"""7 + const.2N^^^-\ (49) 
as long as 7 > 4. If 7 < 4, the first part of the integral is of order 0{N^^/^). 



11 



The second integral of Eq. ([6]), /22, can be estimated in a similar way with the scaling ansatz ^[g},'^] = Qm{NPe). 



The scaling function Gm{x) should have the properties Qm{x) '^^^ and Gm{x) ^-f' const, and should scale as 
N~^e~'^ for e = N~w in or 
gives rise to a contribution 



N foi e = N lo in order to match the neighboring regime, which yields the scaling relation p = -^^ + jq- This 



Jo Jn-^/^o 



'2Nc2e 



4 



_i2._2 logA^ 
= const. 3 ^S'j ^ + const. 4 ^ , (50) 

provided that 77 > 4. If r; < 4 the first part of the integral is 0{N~^) instead. We therefore find the scaling exponent 
of the sample-to-sample-fluctuations with the system size, AF/y ~ N^', to be 



14a 2 /6 \ 3 6 

— — a 



5' 5 2 max(7,4) \5 / ' 10 5max(?7,4) 



(51) 



As decribed in section[3l we have another, slightly different route to chaos and will use it now to check for consitency 
with the above result. From Eq. ([7]) we take the first integral In = m'^ Jq° defi{e){[q'^] — [q^]^) and use the same 
scaling ansatz as above for equation and obtain contributions of the form 

In = r de/i(e)([g4] - [q^f) - const.siV-"-^^?-")-^^ -I- const.eN^^/^ + ... (52) 
Jo 

This integral is positive and has a negative prefactor in Eq. ([7]). The fluctuations can not be negative, however, 
therefore the leading order of this term must be compensated by the second integral /12. We use the scaling function 
GmiN^e) again, which yields 

Ii2=m / degi{e)[q^] const.rN^^^^ + const.giV""^/^". (53) 
Jo 



This second integral, together with its leading prefactor N from Eq. (O, must at least cancel the term of order ©(iV^/^) 
which is contained in N'^In. This is only possible if ?y > |. Hence we obtain a limit on 77 and no contradiction to the 
result derived above 



V. CONCLUSION 



The purpose of this work was to calculate the flnite size scaling of the sample-to-sample fluctuations of the free 
energy in the m-component vector spin glass in the limit of large to. The result is Eq. (j5ip . Although this equation 
looks unpromising at flrst sight, it is in fact very informative. First, we have the solid result /i > ^. Second, from our 
numerical work we know that a > | . But under this condition, the second term in the max function in Eq. (jSip is < ^ 
and can thus simply be omitted. Third, the last term in the max function is < such that we obtain ^ < fJ. < jq- 
Even better, the third term is greater than ^ only for rj > 12, which seems an unlikely large value. We therefore 
conjecture that = ^ is in fact the exact answer. But be that as it may, the exponent 77 could easily be measured in 
simulations, and work along these lines is in progress. 

Since our result is not an exact mathematical proof, we will summarize here the main assumptions on which it 
rests because they may have been obscured by the technicalities. The first ingredient is the connection between the 
fluctuations and bond chaos of the link overlap, Eqs. © and ([7]). They are mathematically exact equalities and pose 
no problem. The second ingredient is the calculation of bond chaos. This is done using large deviation statistics and 
replica theory. We believe that replica theory in principle gives the correct results. It became apparent, however, that 
for the TO component spin glass the "small" deviations play a crucial role. The small deviations statistics are caused 
by finite size corrections of the free energy, which can not be calculated within replica theory. However, by reference 
to earlier resultsSS, we know at least the finite size scaling of the free energy corrections and can thus estimate the 
point where our large deviation calculation becomes valid. The region of the small deviations is then covered by a 
scaling ansatz, which is assumed to cross over smoothly to the region of the large deviations. The introduction of the 
scaling ansatz also introduced a number of unknown exponents. However, only three of these exponents are actually 
relevant for our results, and one of them, a, has been measured experimentally. Moreover, the precise values of the 
exponents are largely irrelevant: if a > | (as we have checked numerically) and 77 < 12, then = ^ is exact. 
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The result /J. = ^ (or even the range < /i < ^) is very interesting because it demonstrates that the large-m model 
is fundamentally different from the spherical model^S^ even though their free energies are identical^^. While the width 
of the distribution of ground state energies of the spherical spin glass scales as N^^^, this behavior is definitely ruled 
out by our results. 

It would be interesting to see whether a result similar to ours could be obtained using the methods of Ref. [111. We 
believe a corresponding replica calculation for the large-m model ought to be feasible. 
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